load Claybank_mc_may112011.txt
d= Claybank_mc_may112011;
%be sure to modify the legend below

x=1:length(d(:,1))  %number of records
plot(d(:,6),'r-')        %turbidity from Turbidometer 1
hold on
plot(d(:,8))   %trubidity from Turbidometer 2
legend('S4912 Core A', 'S4909 Core B')             %Need to identify the cores
[AX,H1,H2] = plotyy(x,d(:,8),x,d(:,7), 'plot');   %disk rotation (RPM) from sensor 2
xlabel('record number', 'fontsize',14)
ylabel('turbidity (NTU)','fontsize',14)
set(get(AX(2),'Ylabel'),'String','disk rotation (RPM)', 'fontsize', 14)
title('YR110511')

% puts the stress values where you want them

gtext({'0.01 PA'},'Fontsize',12)
gtext({'0.05 PA'},'Fontsize',12)
gtext({'0.1 PA'},'Fontsize',12)
gtext({'0.2 PA'},'Fontsize',12)
gtext({'0.3 PA'},'Fontsize',12)
gtext({'0.45 PA'},'Fontsize',12)
gtext({'0.56 PA'},'Fontsize',12)




